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ABSTRACT 


Supersonic combustion Ramjets (Scramjets) are used for h>personic flight speeds where 
the combustion takes place supersonically. For Scramjets, Literature shows that the 
Method of characteristics (MOC) is a good technique for designing a proper nozzle 
contour. A general ‘C’ language computer code is developed for designing a Scramjet 
nozzle contour using the Method of Characteristics numerical technique Initially a 
nozzle contour is designed for parallel flow at the exit. The effect of back pressure on 
nozzle exit at different altitude is evaluated. The variation of exit parameters and nozzle 
contour for different assumed values of constant-gamma and variable-gamma throughout 
the nozzle is studied for air. Variable-gamma is assumed as a function of temperature and 
chemical composition of the exhaust mixture. For different fuels, the species 
compositions of the exhaust mixture are calculated using stoichiometric equations. The 
species concentration is assumed constant throughout the length of the nozzle. The 
gamma variation after each characteristic is evaluated which gives the correct nozzle 
contour for Scramjet nozzle. The developed code can be used for both two-dimensional 
and axi-symmetric nozzles. 

The results obtained from the MOC are verified with the available results in the literature. 
The factors that affect the nozzle performance are analyzed and calculated at different 
axial locations. The final nozzle contour is obtained by truncating the length at particular 
location where the performance parameters are optimum. From analysis, it is found that 
60% of stream thrust function value gives a reasonable length. Force and moment 
coefficients along the nozzle contour at various stations are computed which is used for 
controlling the vehicle where the aerodynamic controls are ineffective. Flow through the 
designed nozzle is analyzed with an available Euler’s code. The exit parameters obtained 
from the Euler’s code are in good agreement with the MOC results. 
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CHAPTER 1 


INTRODUCTION 

Propulsion is a science of designing an engine to propel a vehicle forward or upwards. In 
aviation, propulsion is classified into two types namely, air breathing and rocket 
propulsion. Both the types of propulsions, work on the same principle of pushing high 
velocity exhaust gases out at the back end, but are different in one important way. An air- 
breathing engine uses the atmospheric air as oxidizer, which makes it usable only within 
the atmosphere. A rocket engine is made for traveling in space. It would carry both fuel 
and oxidizer. So the specific impulse of the air-breathing engine is high compared to the 
rocket engine. Turbo jet engine has a specific impulse of around 3000 s, which is 
efficient up to Mach number 2.5. But rocket engines have only around 400 s. Turbo jet is 
useful for a particular case when high acceleration and high speed is not important, but 
most missile applications require high thrust and acceleration. Ramjets are attractive, in 
the high supersonic speed range, but initiation of this speeding up process requires a 
booster. If Mach number is more than 5, Scramjets (Supersonic combustion Ramjet) are 
preferred. 

Ramjets are used in the flight Mach number ranging from 3 to 6. Both Ramjets and 
Scramjets need not be symmetric about the centre line of the engine, because they do not 
have any rotary mechanism. It is often convenient that the bottom surface of the vehicle 
Itself is acting as an internal boundary of the engine. When the Mach number goes 
beyond 5, Scramjets are preferred because the losses in the Ramjets are quite large. 

Scramjets are normally preferred for the flight Mach number range of 5-20. When the 
airflow is decreased by either Scramjet or Ramjet, both the relative velocity and kinetic 
energy of the incoming flow is reduced. This increases the internal energy of the flow 
coming into the^combustor. For Mach numbers above 6, this effect is predominant. It is 
no longer advantageous to decelerate the flow to subsonic speed, as the inlet static 
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pressure increases to very high value This leads to the failure of burner structure, 
excessive wall heat transfer, dissociation of large fraction of available chemical energy, 
and pressure loss in the combustion chamber. So for hypersonic Mach numbers, 
Scramjets are preferred. 

1.1 Nozzles 

The function of the expansion system is to accelerate the flow from the burner exit static 
pressure to local atmospheric pressure over the entire range of vehicle operation in a 
controllable and reliable manner with maximum performance. This is accomplished by 
converting the thermal energy of the hot combustion chamber gases into kinetic energy. 
In simpler terms, nozzle is the component of a rocket or air-breathing engine that 
produces thrust. 

1.2 Types of Nozzles 

There are three major and basic types of nozzles used in aircraft applications namely 
conical nozzles, bell or contour nozzles and annular or plug nozzles. All other nozzles are 
modified from one of these basic types. 

1.2.1 Conical nozzle 

The conical nozzle was often used in early rocket applications because of its simplicity 
and ease of construction. The cone gets its name from the fact that the walls diverge at a 
constant angle. A small angle produces greater thrust, because it maximizes the axial 
component of exit velocity, but minimizing the angle makes the nozzle longer and 
heavier, which leads to complexities while building them. At the other extreme, size and 
weight are minimized by a large nozzle wall angle. Unfortunately, large angles reduce 
performance at low altitude, because the high ambient pressure causes overexpansion and 
flow separation. Here, the exit flow can be achieved with some angularity, which reduces 
the axial thrust of the vehicle. 
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Fig. 1.1 Conical Nozzle 


1.2.2 Bell Nozzle 

This IS the most commonly used nozzle shape, offering significant advantages over the 
conical nozzle, both in size and performance. Near the throat, the nozzle diverges at a 
relatively large angle but the degree of divergence tapers off further downstream of the 
nozzle. Near the nozzle exit, the divergence angle is very small. In this way it minimizes 
weight while maximizing performance. The most important design issue is to contour the 
nozzle to avoid oblique shocks and to maximize performance. However, the final bell 
shape will only be the optimum at a particular altitude. 



Fig. 1.2 Bell Nozzle 

1.2.3 Single expansion ramp nozzle 
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Single expansion ramp nozzles are mainly used for hypersonic air breathing vehicles and 
missiles where the drag is an important performance parameter. In these nozzles, the 
bottom surface of the nozzle is removed because this portion produces large amount of 
drag compare to thrust. The expansion is obtained from simple sharp comer placed at the 
inlet of the nozzle. These nozzles are more efficient as compared to conical nozzles. Fig. 
1 .3 shows one such nozzle. 



Poo 


Fig. 1.3 Single expansion ramp nozzle 


1.2.4 Annular Nozzles 

The annular nozzle, also known, as the plug or "altitude-compensating" nozzle is the least 
employed of those discussed due to its greater complexity. The term "annular" refers to 
the fact that combustion occurs along a ring, or annulus, around the base of the nozzle 
"Plug" refers to the center body that blocks the flow from what would be the center 
portion of a traditional nozzle. "Altitude-compensating" is sometimes used to describe 
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these nozzles since that is their primary advantage. When considering an annular nozzle, 
the area of the center body (A plug) must also be taken into account. 
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Fig. 1.4 Plug Nozzle 


1.3 Application of Nozzles 

Aero spike nozzle is mostly used for rocket applications. This nozzle is also called 
altitude compensation nozzle. Conical nozzle is also used for both rocket and hypersonic 
vehicles. Contour nozzle is mainly used for supersonic wind tunnels to get parallel flow 
in the test section. Recently they are being considered for hypersonic air breathing 
propulsion to increase the thrust. In order to design such a nozzle, initially a hypersonic 
parallel flow nozzle is designed using Method of characteristics numerical technique, and 
then this nozzle is truncated for rocket application, using various performance 
parameters. It is important to note that all conventional supersonic aircrafts have contour 
nozzle. 

1.4 Nozzle Design 

The function of the ideal nozzle is to provide parallel flow at the nozzle exit for any 
nozzle exit conditions. Any divergence in the nozzle leads to loss in its performance. 
Hence the contour nozzles are preferred for aircraft and wind tunnel applications In 
contour nozzle the exit uniform flow depends on exit Mach angle. If exit Mach number is 
greater, then the length of the nozzle also increases to very high value. In hypersonic 
vehicle the length of the nozzle is an important parameter. If the length increases then the 
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drag produced by this length itself is very large, so there will be more drag in addition to 
parallel flow and high thrust. The reduction in length is accomplished by using sharp 
comer expansion instead of continuous expansion. Method of characteristics is widely 
used for more practical nozzle design. This method is also used to design a nozzle for 
thennally perfect and chemically reacting flow with both equilibrium and non- 
equihbnum effects. This is a simple numerical method for designing a nozzle contour. In 
this method the flow is inviscid and irrotational. 

Engineers in the conceptual air breathing vehicle design need to be able to assess engine 
performance at overy point in a vehicle’s trajectory. Unlike rocket engines, where 
propulsion system is dependent only on the vehicles altitude, air breathing propulsion 
system performance is strongly dependent on the flight path (velocity, dynamic pressure) 
of the vehicle as well. 

In the rocket or air breathing missiles, the nozzle not only produces the thrust, but also 
produces a normal force on the nozzle surface of vehicle. This is an additional source of 
lift available on the vehicle. The pressure distribution along the nozzle surface greatly 
affects the pitching moment of the vehicle. This pitching moment is produced because of 
the expansion surface of the nozzle. These moments influence the stability of the vehicle. 
The force and the moment coefficients should be calculated precisely so that the 
equilibrium flight is sustained when the engine power is switched off These moments are 
also used for adjusting the path of the missile. This phenonienon is called thrust 
vectoring. Here the hypersonic nozzle contour is designed and optimized for single 
expansion ramp nozzle, which is used in air breathing missile propulsion. Force and 
moment coefficient and the overall thrust of the vehicle are also calculated. The length of 
the contour nozzle is also reduced based on different parameters. Finally the flow 
angularity also calculated. Details of the analysis are given in the following chapters. 
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CHAPTER 2 


LITERATURE SURVEY 

Encountering propulsion problems are complex m real life, as interplay of heat and mass 
transfer, fluid mechanics, chemical kinetics, thermodynamics and turbulence occurs. 
Such problems are not tractable analytically, as the non-linear nature of governing partial 
differential equations involved in the flow makes it very difficult to get closed form 
solution. This leads to the existence of numerical methods. Method of characteristics is 
one such numerical technique, which is applicable for both super sonic and hypersonic 
nozzle design. 

With an increasing demand for the reduced cost of space and missile activities, 
researchers are forced to investigate new technologies, such as increasing fuel economy. 
Greater fuel efficiency can be achieved in air breathing missiles and spacecrafts with the 
design of better propulsion and suitable selection of nozzle contours. In order to achieve 
maximum thrust of the propulsion system, recent researches focus on proper designing of 
the after body (nozzle). A push for increased technology in this field has lead to the 
development of optimum air breathing vehicle nozzle contours that can operate 
effectively and reliably at high altitudes. 

2.1 Existing methods for Nozzle design 

In the early years of rocket mid missile nozzle design, conical nozzles were designed 
because methodology for designing efficient contour nozzles did not exist. Although they 
were necessary at that time, conical nozzles were an inadequate solution. Such nozzles 
resulted in performance losses, because they did not produce axially directed exhaust. 
The angularity in the conical nozzle leads to incorrect nozzle design for wind tunnel 
applications. In some cases, these nozzles were excessively heavy because of their length, 
which leads to enormous drag. This is one of the main reasons for hypersonic flights were 
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not possible in early days. Therefore, the design of efficient nozzle and parallel flow 
exhaust was a challenging task. One of the major problems in designing a supersonic 
wind tunnel is the determination of appropnate nozzle contour. Eradicating this problem 
would yield a parallel and uniform flow in the test section. 

Short and efficient contour nozzles were finally developed when the truncated perfect 
nozzle approach was pioneered. In this approach, a wind tunnel nozzle was designed for 
uniform flow and then truncated at a much shorter length. Using this approach, Ahlberg 
et al [1] presented some results that were obtained for supersonic nozzle. Hoffman [2] 
discussed some variation of the original method. With in the same penod, an alternate 
approach to this problem, based on the calculus of variations, was formulated. This 
method was improved by Rao [3]. According to the procedure outlined by Rao, the 
length, ambient pressure, and flow properties in the immediate vicinity of the throat, are 
the governing conditions under which thrust is maximized. 

This classical process assumes that an isentropic flow exists in the nozzle. The nozzle 
contour is constructed using Method of characteristics to give the desired flow al the exit. 
Historically, Rao’s method has been considered the best, as the other classical approaches 
fail to surpass its performance. However, the drawback in Rao’s method is its 
dependence on an inviscid and calorifically perfect gas. Classical procedures are based on 
an approximation of viscous effects. These methods rely on an inviscid design. 
Therefore, this method is inapplicable for high temperature gas flow in the hypersonic 
nozzle design and certain correction must be made for types of fuels used in the aircraft 
propulsion devices. 

Formulae required for calculating the species specific heat capacity and specific heat ratio 
at very high temperature were proposed by Oussama Rizkalla [4]. Computational studies 
on calculated chemical and vibrational Non-equilibrium Effects in hypersonic flows were 
also done. In order to study vibrational effects of both equilibrium and Non-equilibrium 
conditions in hypersonic flow, some results from quantum mechanics were used. 
Calculations for both force and moment coefficients for axial and normal directions were 
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also made for a Mach number 20 flight condition. From these analyses, it was concluded 
that the vibrational relaxation time appear to be quite short at these conditions. Hence, the 
flow appeared to be in vibrational equilibrium throughout its passage through the nozzle. 

The vehicle trim requirements were also determined from calculations that include both 
finite rate chemistry and viscous effects. Beyond a certain axial location, which is 
considerably near the upstream of the nozzle exit station, certain amount of additional 
axial thrust was obtained. This additional thrust is because of the added expansion surface 
m the after body of the nozzle. However, the added expansion surface significantly 
affects the normal force coefficient and its associated pitching moment. The pitching 
moment will significantly affect the vehicle tnm requirements. Therefore, care should be 
taken while calculating the normal force and pitching moment. 

A simple method for calculating the properties of combustion products was given by 
Guha [5]. The method is generic: the same equations apply to any fuel for a specific 
temperature range. This offers best accuracy comparable with the most elaborate method 
available but at a significantly reduced computational cost. The numerical methods 
presented here are for lean mixtures (excess air). For nch mixtures, when not all the fuel 
is burned, the chemistry of the combustion process can be complex. However, if the 
composition of the products of combustion is known, then this simple methodology can 
be extended to rich mixtures. 

The MOC method would offer considerable savings in computational time with no 
significant loss of accuracy. So, this method is used to calculate th? specific heat capacity 
of any combustion species once the temperature in the combustion chamber is known. In 
the current nozzle design, it is assumed that complete combustion takes place in the 
combustion chamber. The maximum combustion temperature is attained when 
hydrocarbon fuels are mixed with just enough air so that complete conversion of the 
hydrocarbons to water vapor and carbon dioxide occurs. 
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“An Introduction to Numerical Techniques for Nonlinear Supersonic Flow” was given by 
John D. Anderson [6], where both Method of characteristics and finite difference 
numerical methods were discussed. Meanwhile, newer methods to solve the nonlinear 
partial differential equations that exist in the potential flow were given by Michel Saad 
[7] These nonlinear partial differential equations were solved using linearized potential 
flow (small perturbation theory) and Method of characteristics technique. A detailed 
procedure of Method of characteristics design was also discussed. 

A simple procedure for designing a low supersonic contour nozzle design was described 
in “Gas Dynamics” by E. RathaKrishnan [8] where the gas is assumed calonfically 
perfect. However, this cannot be used for supersonic and hypersonic nozzle design. Here, 
the present design is for thermally perfect gas where the chemical composition of the 
exhaust was frozen. 

2.2 Scope of the work 

Method of characteristics numerical method is used here for designing a hypersonic 
Single expansion ramp nozzle. The minimum length hypersonic contour nozzle for wind 
tunnel application is designed and the length of the nozzle is optimized for missile 
application using “C” language computer code. The nozzle contour is designed for 
different Mach number and altitude conditions. The length of the nozzle is minimized 
using different design parameters such as stream thrust function, force and moment 
coefficients. 

Here, the combustion is assumed to be complete. So, using stoichiometric equation the 
exact species composition is found. From the given inlet, temperature in the nozzle, 
specific heat capacity and specific heat ratio are foimd using Guhd’s generic method [5]. 
The chemical composition of the nozzle exit is constant and the temperatures at different 
axial location are evaluated from simple isentropic relations. After this, nozzle contour is 
designed using Method of characteristics numerical method. Provisions are also provided 
for different altitude atmospheric condition calculations. These formulae are taken from 
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NASA Online tutorial [27]. This design is correctly expanded nozzle design. Force and 
moment coefficients are also calculated along the length of the nozzle using the formulae 
given by Oussama Rizkalla [4]. 


Stream thrust function from the nozzle inlet to exit is calculated. Stream thrust function is 
a non-dimensional function that is used to find the mass flow rate specific thrust of the ' 
vehicle. Momentum and pressure thrust along the nozzle is also calculated. Axial and 
normal force and moment coefficients are computed for entire length of the nozzle. The 
length of the hypersonic nozzle is optimized from the above-mentioned parameters. The 
normal force coefficient leads to additional lift on the nozzle section, which produces the 
pitching moment. The knowledge of pitching moment is useful for controlling the vehicle 
where the aerodynamic controls are ineffective. 

The results obtained from this design is verified with Ref [4] available data and verified 
with Euler code where the flow is inviscid. In Euler code, the nozzle geometry is 
designed using, the results obtained from present work and the results obtained from 
Euler code are compared with MOC results. 
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CHAPTER 3 


THEORY 


Analytical solutions of the equations describing fluid flow have been found only for 
simplified problems. No general analytical solution has been found yet for non linear 
partial differential equations. For this reason numencal techniques have been used for 
solving non linear problems. One such technique is the Method of characteristics (MOC), 
used for ordinary two dimensional exhaust nozzles, which form the basis for design of 
Ramjet and Scramjet exhaust systems. 

The MOC method can be applied to one dimensional unsteady flow, two and three 
dimensional steady flows and axi-symmetric steady flow .In all the cases flow must be 
irrotational and isentropic. The computer code developed here can be used for both two 
dimensional and axi-symmetric contour nozzles design. 

The Method of characteristics indicates the properties of a perfect gas m which 
continuous waves of small but finite amplitude are present as the gas flow supersonically, 
irrotationally, and shock free in any passage. 

From the equation of law of conservation of mass for unsteady flow, 

dp ^djpu) ^d{pv) ^d{pw) 

dt dx dy dz 


For steady (independent of time) flow, we can neglect the first term, 


du dv dw 
dx dy dz 


( 


dp dp dp 

y_£_ + v—t- + ^ 

dx dy dz 


=0 


( 1 . 2 ) 
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These components can be expressed m ternis of the velocity potential as: 


u 


dx 




d(j) 

dy 

^_d(p _ 
dz 




The continuity equation in terms of velocity potential may then be written as: 


K + + - 

p 


^ , dp , dp dp^ 
(b — + 6 — + d> — 
OX ay dz 


=0 


(1.3) 


From the momentum equation, assuming that only pressure forces are present (neglecting 
Body force and viscous force). 


— + 
P 


dV^ 

2 


= 0 


But = U^ + V^ + ^ ^ .V + y + (j) ^ z 

Substituting this value in the above equation, we get 


(1.4) 


— + -[^\+(l)\+f-z)=0 (1.5) 

P 2 

We know that, for isentropic flow: 

dp 

Velocity of sound (c) is applicable for any flow (thermally perfect and chemically 
reacting flow). 

Eliminating d p from equation (5): 



Partial differentiation of this equation with respect to the three coordinate axes gives: 


d Xf c 

^ = -4 + («,(!>,,+ A «>„) 

ay c" 

^ = -^ + +«>.(*=) 


Substituting these values in equation (3), 


<l>n + <l>,y + - 



+ <l>y{^x<f>xy + <l>y<l>yy+(l>2<l>yz) 

^<f>z^x^xx + <t>yAy + ) 


= 0 


+ <l>yy + <l>: 


“T [^^’‘^xx + <l>^y<l>yy + ) 

c 

-\k>x<t>y<^xy + (l>y<t>2<l>yz )= 0 

C 


(1.7) 


This is the differential form of potential function in three dimensions. 

For two dimensional flow, 

«»„+(»„. -4(«s“<«>»+f»^(»„)-4{«A)=o 0.8) 

C C 


Rearranging this equation, we get 



This equation is similar to second order partial differential equation with two independent 
variables. 
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( 1 . 10 ) 




dx- 
Where, 


dxy dy^ 

a=,-4 




The terms A, B and C are functions of both the independent variables x, y and the 
dependent vanables u, v. An important point to be noted here is that the velocity of sound 
c is not a vanable. This can be expressed in terms of dependent variables. 


There are many boundary value problems which involve partial differential equations 
Only few of these equations can be solved by analytical methods. In most of the cases we 
depend on the numerical solution of such partial differential equations. Of the various 
numerical methods available for solving these equations, the method of finite difference 
is commonly used. Computational fluid dynamics is one of such method used for 
aerospace applications. 


Computational Fluid Dynamics (CFD) is a new third dimension in aerodynamics, 
complementing the previous dimensions of both pure experimental and pure theory. 
Basically, CFD is the art of replacing the governing partial differential equations of fluid 
flow with numbers, and advancing these numbers in space and time to obtain a final 
numerical description of the complete flow field of interest. 

In this method, the derivatives appearing in the equation and the boundary conditions are 
replaced by their finite difference approximations. Then the given equation is changed to 
a difference equation which is solved by iterative procedure. This is slow but gives good 
results of boundary value problems. An added advantage of this method is that the 
computation can be done by electronic computer. 
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MOC IS simple and efficient method for inviscid supersonic and hypersonic nozzle 
design over computational fluid dynamics. This method can be applied to one 
dimensional unsteady flow, two and three dimensional steady flows and axi-symmetric 
steady flow. 

Equation (10) has three possible solutions, 

(i) - AC <0 Elliptic equation. 

This equation has two imaginary roots. There is no specific direction for 
propagation of waves, so the wave will propagate like a spherical wave. 

(li) - AC =0 Parabolic equation. 

This equation has two identical real roots. The wave will propagate in only one 
direction. 


(iii) B^ - AC >0 Parabolic equation. 

This equation has two real and distinct roots. The wave will propagate in two 
distinct directions. 

Now consider the equation (10), 

2 V- 
B^-AC=^-1 
c“ 

For Subsonic flow V<c, then B^-AC <0, which is nothing but the elliptic equation. 

For Sonic flow V=c, then B^-AC =0, which is parabolic equation. 

For Supersonic flow V>c then B^-AC >0, which is hyperbolic equation. 

The problem under consideration is Hypersonic nozzle design, where V>c. 

Hence there are two distinct real roots. These solutions are associated with two families 
of characteristic curves in the x-y plane, which are known as physical characteristics. 
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As an example, consider a uniform isentropic supersonic flow past a convex comer 
shown below, 



Fig. 3.1 Flow over convex corner 

The Mach waves are generated at the comer, so that the gas flows uniformly, but at 
different velocity, while crossing each wave. Since the flow is shock free, the velocity is 

continuous, so that the first derivatives of the velocity potential — and — are 

dx dy 

continuous The gas expands as it flows around the comer, therefore increasing the 
velocity when the gas reaches the region downstream of the wave. This means that there 
is a discontinuity in the velocity gradient at the Mach wave. 

Since the velocity components are continuous, the changes in velocity components are 
given as: 
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( 1 . 11 ) 


, du , 8u , 
du = — dx + — ay 
dx dy 

and 

dv dv 

dv = — dx H ay 

dx dy 

However, the condition of irrotationality implies that = 0 

dy dx 


From this 


, du , dv . 
dv = — dx + — ay 
dx dy 


Equation (10) can be written as, 
A 

Rearrange the Equation (11), 
dx 
and 


(du') 

+ 2B 


+ C 


— 

[dxj 






= 0 


f du'^ 

+ d y 

(du') 

+ 0 







= du 


0 


du 

dx 


+ dy 


y^yj 


+ dy 


^dv^ 


\^yj 


= dv 


This set of equation can be solved for each of the three unknown. 
These equations are solved using Cramer’s rule, 


d~(/> 

dxdy 


'-4 



0 




c" 

dx 

du 

0 

0 

dv 

dy 

i-!d 

7 

c 

-2uv 

>4 

c 


dx 

dy 

0 

0 

dx 

dy 


D 


(1 12 ) 
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Consider the fig below. 



Fig. 3.2 Flow description 


The derivative 


d-(p 

dxdy 


has a specific value at point A. Equation (12) gives the value 


solution for 


d'(j) 

dxdy 


for an arbitrary choice of dx and dy. The combination of dx and 


dy define arbitrary direction ds away from point A. This direction is different from the 
stream line direction going through point A. The differential du and dv represent the 
change in velocity that takes place due to increments in dx and dy. Hence, although the 
choice of dx and dy is arbitrary, the value of du and dv must correspond to this choice 
No matter the value of dx and dy are arbitrarily chosen, the corresponding value of du 


and dv will always ensure obtaining the same value of 


ay 

dxdy 


at point A. 


The single exception to the above comments occurs when dx and dy are chosen so that 

aV 


D=0 in equation (12). In this case 


dxdy 


is not defined. This situation will occur for a 


specific direction ds away from point A, defined for that specific combination of dx and 
dy for which D=0. 
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ay 

dxdy D 0 
d-(j> 

Here, — — is an indeterminate form. The important conclusion here is that there is 

dxdy 

d^ d) 

some direction (or directions) through point A along which — is indetemiinate. This 

dxdy 

implies that the derivatives of the flow variables are indeterminate along these lines. 
Hence, we have proved that lines do exist in the flow field along which derivatives of the 
flow variables are indeterminate .These lines are called characteristics lines. 

Setting the numerator equal to zero, gives: 


A (d yf - 2B dx dy + C (dxf =0 

{dx)^ =0 


f 2 ^ 

{dyY + '^^^dxdyy 

f 2 \ 

1-4 

V ^ J 




f 2 ^ 

(dy'] 

2 2uv 


■f 

f ^ 

1_2_ 

1 ■ c J 

^dx j 

, 

char 

\dx ^ 

char 



(1.13) 






(1.14) 



(1.15) 


From above figure, the velocity components u and v are: 
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u = Vcos0 
V = Fsin^ 

Substituting these values in equation (14), we get 


= tan(6'+/i) ; (1.16) 

char 

It states that two characteristics lines run through point A, one line with the slope equal to 
tan(6'-//)and other with the slope equal to tan(^+//). These lines are shovra below. 




Fig. 3.3 Inclination of characteristic curves 

The characteristic lines through point A are simply the right and left running Mach 
waves. Therefore, Each and every point the flow field has two characteristic lines, up and 
down running characteristics. It is very difficult to consider infinite number of 
charactenstics in flow field analysis. For practical applications we deal with a finite 
number of waves. Note that the characteristic lines are curved in space because, 

(1) The local Mach angle depends on the local Mach number, which is a function of x 
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and y 

(2) The local free stream direction 6 varies throughout the flow. 
Now considenng N=0, then 


1-^ 



du dy + 

1-^ 



1 


dv — 1 


(dy^ 

du (l “ v~/c“ ) 

\^dx] 


Substituting the value of 


dx ^ 


char 


from equation (14), we get 



Substituting the value of u and v in term of resultant velocity (V), and after some 
algebraic manipulation, we get: 

de = + 4M^ - 1 ^ 

From energy equation. 


1 = Const3in.t 

v-l 2 


But c = V/M, and there fore: 

; H = Cons tan/ 

2 


22 



Differentiation of this equation with respect to V leads to an expression for velocity 
change as a function of Mach number: 


dV 

V 


dM- 


2M- + 


Now, substituting this value in 


equation (18) and after integration: 


v{m)+ 9 = const 
This can be written as, 

u {M ) + 9 = const = K_ (Along the down running characteristic) 

u{m)- 9 = const =K^ (Along the up running characteristic) 

Note that our compatibility relations are now reduced to algebraic equations which hold 
only along the characteristic lines. In a general, steady inviscid supersonic flow 
compatibility equations are ordinary differential equations. 


These results are used to determine the supersonic flow inside the nozzle and the 
determination of proper wall contour, so that shock wave does not appears inside the 
nozzle. The detailed design procedure is given in the following chapter. 


23 



CHAPTER 4 


METHODOLOGY 

Method of characteristics is a simple and efficient method for supersonic nozzle design. 
This method can be applied to dne-dimensional unsteady flow, two and three- 
dimensional steady flows and axi-symmetnc steady flow. The method of characteristics 
indicates the properties of a perfect gas in which continuous waves of small but finite 
amplitude are present as the gas flows supersonically, irrotationally and shock free in any 
passage. 

The compatibility relation is given by, 

v{m)+ 9 = const =K_ (Along the down running characteristic) ...( 4 . 1 ) 

v{m)- 9 = const = (Along the up running characteristic) .... ( 4 . 2 ) 

Consider the following flow field, the flow properties at different locations are calculated 
using the following procedure. 



Fig. 4.1 Grid points in the flow field 
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With the given geometric coordinates and flow properties along a starting line, a number 
of initial points, such as points a, b, c, etc., are selected. From the values of 0 and y at 
each of these points, numerical value of characteristics Ci and Cn are determined. At point 
d, characteristic Ci passing through point “a” intersects with characteristic Cn passing 
through “b”. The value of 6 and y at point “d” is calculated using the above compatibility 
equations. 

For example, the two characteristic curves, which locate point “d“, are 

(^1 = +^a=^d + 

and 

^ U ~ ~ ~ ~ 

The resulting flow direction 0 and the Prandtl-Mayer function o at point “d” are 
determined by these characteristic values, Ci and Cn 

g k ^ c,-c„ 


and 

'^^b)+{^a ~^b) _ Ci + C„ 

2 ~ 2 

The Mach number and Mach angle corresponding to can be determined. Since the 

flow is isentropic, temperature and pressure at point “d” can be calculated from isentropic 
relations. Similarly, characteristic Ci of point “b” and characteristic Cn of point “c” are 
used to locate point “e”. Point “f ’ has the same Ci characteristic as point “a” and “d”, and 
the same Cn characteristics as points “c” and “e”. The procedure is continued to other 
points in the flow field, and the complete net of characteristic curves describing the flow 
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pattern can be determined. Knowing the inclination, the characteristic lines allow the 
construction of the characteristic net in physical plane. This characteristic net gives the 
physical meaning of the flow. 

4.1 Formulae 

4.1.1 Atmospheric condition calculation 

Nozzle exit Mach number and exit area depend on back pressure (or ambient pressure) of 
the nozzle. This back pressure value is not constant for rocket and missiles Similarly, 
temperature and density also vary with altitude. The variation of static temperature with 
geometric altitude consists of series of straight-line segments either isothermal or 
gradient. The static pressure decreases exponentially with increasing geometnc altitude. 

The variation of these thermodynamic parameters for different altitudes from Ref [27] is 
given below: 


Altitude <11000 m 


Temp 

Pres 


Dens 


= 15.04 -(0.00649*//) 


= 101.29* 


^ temp + 273.1 
^ 288.08 ^ 


Pres 

(0.2869 (Temp + 273. l)) 


(4.3) 


11000 m< Altitude < 25000 m 


Temp = - 56.46 
Pres = 22.65 


Dens 


Pres 

(0.2869 {Temp + 273.1)) 


(4.4) 
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Altitude > 25000 m 


Temp 
Pre^ = 2 488* 


Dens = 


-131.21 + 0.00299*// 

Temp + 273.1 
216.6 
Pre5 

(0.2869 {Temp + 273. l)) 


(4.5) 


4.1.2 Specific heat ratio calculation 

Normally, Specific heat ratio y is assumed constant for many applications, such as design 
of wind tunnel nozzle. In ordinary supersonic wind tunnel, the maximum temperature is 
around 300K. A gas is considered calorifically perfect up to 800 K. If the temperature is 
above 800 K, the vibrational degrees of freedom present in the gas will increase. 
Therefore, if the temperature is greater than 800 K but less than 1750 K, the gas can be 
considered as thermally perfect, which implies that, the gas properties are function of 
temperature alone. When the temperature exceeds 1750 K, the gas properties are 
functions of both temperature and pressure. These kinds of flows are called chemically 
reacting flows, because at this temperature the gas begins to dissociate. O 2 molecules 
dissociate to oxygen atoms and N 2 molecules dissociate to Nitrogen atoms. Due to this 
phenomenon, the amount of heat required to increase the temperature of one Kg mass 
will reduce. Specific heat ratio value ( 7 ) depends on the number of degrees of freedom 
present m the gas. When the temperature exceeds 8000 K, ionization of gas atoms occurs. 
The atoms liberate electrons and this leads to ionized state of the gas. This kind of flow is 
called plasma dynamics. 

In Rocket, the temperature at nozzle inlet is very high. Once the gas expands in the 
nozzle, temperature is reduced to very less value along the length of the nozzle and the 
ratio of static to total temperature margin is reduced to very less value at the exit. Due to 
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this drastic variation in temperature, it is very difficult to use constant y throughout the 
length of the nozzle. 

The mixture specific heat is given by 


Cp=E “.c, 


(4.6) 


Specific heat ratio is given by 





(4.7) 


The maximum combustion temperature occurs when hydrocarbon fuels are mixed with 
just enough air so that all of hydrogen atoms form water vapor (H 2 O) and all the carbon 
atoms form carbon dioxide. This particular nature of fuel and air is represented by a 
general chemical equation for complete combustion called stoichiometric equation. 


This is a given by 


,, + 1 


( 

X + — 

02 

'79 1 

+ 

— ^ X CO 2 "t” ~ H 2 0 + 

79 

f 

x + — 

1 4 J 

\ 

21 ; 

' 2 ' 

21 ' 

1 4j 


iV. 


(4.8) 


m = S m, 


m 


CO. 


= 44x 


m„0 


m. 


= ^18 
2 

= 11 ^ 
~ 21 


y 


X + 

V 


28 


28 



Total mass = 44 x + 18 


2j 21 


X 4- — I 28 


( 4 . 9 ) 


a, 


m, 


m 


( 4 . 10 ) 


C,„„, =- 0.0849 + 0.69384 -0.09326 (^J+ 5.49954k lO"* 

C,„^„ = 7 9472 - 10.19667 + 4.59728 - 0.20549 {t/^J 

C„_ =1 395-18,3139(i:{ggy'’ +38.31l(^gg)-’ -29.3(r{gg)-’ 


This formula can at least be used in the range 200 - 2000 K. 
The specific heat capacity of the mixture 


C„ 


= c. 


^(h 


a 


CO, 


+ Cr 


niio 


a 


H,0 


+ C 


Pn 2 ^^2 


( 4 . 11 ) 


The characteristic gas constant value is a function of species concentration of the mixture. 
This is also calculated in the similar way. 



~ ^COi ■*' ^N2 (4-12) 

The free stream specific heat capacity is also calculated because, this is used as a non- 
dimensional parameter. Here it is assumed that the atmospheric air contains 79 % of 
nitrogen and 21 % of oxygen. 
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Calculation of Cp 

C,,,^ = 1 16975 + 6.28187 x lO'* ’ 


5.5803 



+ 7.4025 



iVf)^ =32x21 


/72yv, =28x79 


m — itIq^ + m 



So, the specific heat ratio (y) at every characteristic zone is calculated from the local 
temperature of the flow. 

4.1.3 Calculation of exit Mach number and exit area 

The exit area and exit Mach number are calculated using isentropic relation. The nozzle 
inlet area, inlet static pressure, inlet static temperature and nozzle inlet Mach number are 
the input parameters of the program. The specific heat ratio value at the nozzle inlet is 
calculated from the inlet temperature. From the law conservation of mass, mass flow rate 
IS constant through out the length of the nozzle. 

Mass flow rate = pA V 

This can be modified using isentropic relations. 


Mass flow rate (m) 


'RT. 


M A 


1 + 




^2(r-i) 


(4.13) 
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From the given pressure, temperature and Mach number we can , calculate the total 
temperature and total pressure at the inlet using isentropic relations. 


T 

^0/ 

T 

OD 



Po; 



r-1 

2 



J 


(4.14) 


(4.15) 


For Isentropic Flow, 

For coirectly expanded flow Pg = ambient pressure. Ambient pressure is calculated from 
the formulae given for different altitudes. The given nozzle design is shock free correctly 
expanded flow nozzle. 


Assuming a value for Gamma, for example =1.124 


P. 


Of 

p.. 




n.-i 




1 + 

V 2 










Pres/ 


Pi 


= 


r,. -I 




r-y 


Pres 


P, 


rc -1 


(4.16) 


The exit temperature ratio is calculated from exit Mach number. Total temperature 
remains constant throughout the length. Hence, the exit temperature is calculated using 
this formula 
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From this, the local specific heat ratio value of the mixture can be calculated 
Using iterative procedure, the exact value of gamma (i;)can be calculated. Once the value 
of gamma is calculated, then velocity of sound value can be calculated. 

Velocity of sound (a) = ^ y RT (4.17) 


Once the velocity of sound at the exit is calculated, the exit velocity can be calculated 
from this formula 


Exit Velocity (F^) = 


Mach angle is calculated from the given formula, 

^ 1 ^ 


Mach angle (ju) = Sin' 




Area ratio at any section in the nozzle is given by the formula, 





Y + 1 

2(r- 0 

ii. - 


f 2 ] 

( 

A. 


{r + ij 
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1 + 




.. ht L 


J_ 

M 


(4.18) 


From known values of M ^ and Ai, the throat area is calculated. Once the throat area is 
calculated using the area ratio relation, the exit area required for calculated exit Mach 
number can be determined. 

4.1.4 Calculation of total deflection angle 


When supersonic flow is turned through a convex comer, the flow expands lesulting in 
an increase in velocity and drop in pressure .This expansion takes an angle up to 130.5 
degree. The change in properties cannot occur abmptly across an expansion fan but rather 
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gradually across a series of waves emanating at the comer. Therefore, at the comer, flow 
IS not isentropic. If the surface is a continuous expansion surface then flow is isentropic. 

The Prandtl-Mayer function is given by the formula, 

“ tan-'^^(Af=-l)-tan-' V(K’- i) (4.19) 

Prandtl-Mayer function is a similarity parameter. The derivation of this function is given 
in ref [7]. This is based on the assumption that the flow-turning angle is small. 

Knowing the value of inlet and exit Mach number, the total flow-turning angle can be 
calculated. The design is a minimum length exhaust nozzle. This is achieved by placing 
sharp comers that generates centered simple or Prandtl - Mayer expansion fans at the 
nozzle entry, and also achieved by absorbing or canceling each of the charactenstics upon 
their arrival at the opposite nozzle boundary, by means of change in slope of equal 
magnitude and sign. 

Total flow turning from Nozzle inlet to exit is given by. 


o 



tan"' -y/ (m/ - 1) - tan"' ^ (m,' -l) 


....(4.20) 


The amount of flow turning caused by the initial expansion fan must be precisely equal to 
the amount of opposite flow turning along the reminder of the boundary. Therefore, each 
comer must be equal to half of the total flow-turning angle associated with the expansion 
process. Drop in static pressure is one of the parameters to decide the length of the 
nozzle. This parameter is a function of exit Mach number. 

Static Pressure ratio is given by. 
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p. 


1 + 

2 


r-\ 


( 4 . 21 ) 


The ratio of Height at any given location is given by, 

r + 1 

M. 


IL 

H. 




1 + 


2(y-l) 


( 4 . 22 ) 


This is denved from simple isentropic relation. 


Stream thrust function is one of the parameter .to measure the total expansion of the flow. 
The stream thrust function is the parameter, which leads to the determination of mass 
flow rate specific thrust. This is important parameter for performance evaluation. 


The stream thrust function is given by, 
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1 + y M/ 
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4.1.5 Calculation of x and y co-ordinates 





(a^ 




1 

^A, j 

l-\ ^ 


tan(mrn, ^ ) 


( 4 . 24 ) 



( 4 . 25 ) 
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Unifomi flow region length is given by. 



In contour nozzle, the flow at the exit is unifomi. Therefore, after the last characteristic 
the flow must be uniform. The uniform flow length depends on last characteristic Mach 
number. Therefore ninety percent of the length of the nozzle is last charactenstic length 
in hypersonic speeds. Hence, we can truncate the bottom portion of the nozzle after the 
last characteristic wave. The required flap length is given by. 

Flap length = x - 

Still this length is too large for hypersonic vehicle. The length of the nozzle is finally 
optimized according to the values of drop in static pressure and gain in stream thmst 
function. 

4.1.6 Intermediate point calculation 

The intermediate points (characteristic intersection points) of characteristics are 
calculated using simple analytical geometry. From this the idea of charactenstic and the 
physical appearance of the characteristics are known 
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For second down running characteristic, the first intermediate point (1, 2) is calculated 
from the following method. 

From two-pomt fonnula, second down running characteristic equation is written by, 

y2-yi 

xi = 0 and yi = h_i 

-xh_i =X 2 {y-h_i) 

xh _i+X 2 y — X 2 h _i =0 (4-27) 

From one point and slope fonnula, first up running characteristic equation is written by, 

y-y, =mx-mx^ 

ffix-y+y, -/MX, =0 (4.28) 

These two simultaneous equations are solved using Cramer’s rule. The intermediate 
points at remaining locations are also calculated from this formula. Here, the assumption 
is that the characteristic is straight line in between two points. 

4.1.7 Thrust calculation 

The main objective of the nozzle design is to increase the thrust from the available 
pressure at the nozzle inlet of the vehicle. Thrust of the vehicle has two terms. (1) 
Momentum thrust is the main thrust in hypersonic vehicle. The main aim here is to 
increase the momentum as large as possible because, at the exit lot of available energy is 
wasted in the form pressure and temperature. This can be achieved by proper design of 
the nozzle contour. (2) Pressure thrust is an additional thrust in the nozzle in the case of 
under expanded nozzle. 
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Pressure thrust value is zero for perfectly designed nozzle. This means, all the available 
pressure energy is converted in to kinetic energy but this is not possible in hypersonic 
vehicle because the length required for correctly expanded flow is too large The length 
of the nozzle depends on various parameters. Therefore, the pressure thrust is also an 
important parameter in hypersonic vehicle unlike subsonic vehicle. 


In general, thrust of the engine is given by 


T=| k +'’',)K-rx. K |+A(n -P.) 


■> ◄- 


(4.29) 


Momentum thrust 


Pressure Thrust 


Static thrust of the engine is given by, 

T= (m„ +m^)V^ +A{Pe-Pa) 


Assume fuel air ratio is small, 

T=m, V,+A,{P,-P.) 


(4.30) 


4.1.8 Force coefficient calculation 

The extended surface in the nozzle exit produces axial and normal force at the exit of the 
nozzle. This is because of the flow angularity in the flow in the nozzle. These forces 
produce axial and normal moments in the nozzle exit. In Rocket and missile, these 
moments are very important to control the vehicle where aerodynamic controls are 
ineffective. In satellite, small rockets are used to put the satellites in correct orbit where 
the same principle is used. So, the value of these coefficients must be calculated precisely 
because this will lead to unstable vehicle. 
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Axial Force Coefficient is given by, 


=[■ {P-Pj)^tz.n0^,dx (4.31) 

Normal Force Coefficient along the wall is given by the formula, 

= [■' (P-P„) dx (4.32) 

Normal Force Coefficient along the cowl (Flap) is given by the formula, 



Axial Moment Coefficient is given by, 

C^. = {P -Pj^ tme,{y-l)dx (4.34) 

Normal Moment Coefficient along the wall is given by the formula, 

C,, = r (P-P„) xdx (4-35) 

y Hall Jb ^ W 

Normal Moment Coefficient along the cowl is given by the formula, 

C,, = r {P-PS) (4.36) 

Net Normal Moment Coefficient is then calculated from this formula. 
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( 4 . 37 ) 


Here, pressure is normalized using dynamic pressure at exit. 

Length (x) is normalized using throat area 

Finally, force and moment coefficients are calculated for a given configuration. From 
analysis, it was found that this value is large at the initial portion of the nozzle after that 
this value is small because the static pressure of the nozzle is decreased to local ambient 
condition pressure. That is discussed in the next chapter. 
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CHAPTER 5 


RESULTS AND DISCUSSION 

5.1 Validation of results 

The results obtained from Method of characteristics (MOC) numencal solution is verified 
with available h 3 /personic data given in Ref [4]. The results presented in this section are 
for two-dimensional hypersonic nozzles. This code is also used for axi-symmetnc flow. 
This code incorporates the option to use different fuels in combustion chamber. Here it is 
assumed that the combustion-taking place in the combustor is complete combustion. 

In Ref [4], the nozzle contour is designed using MOC and the flow is analyzed using 
Euler’s code. In addition, the flow is considered as chemically reacting with vibrational 
non-equilibrium. 

Input parameters for compared results are. 

Nozzle inlet pressure = 154100 N/m^. 

Nozzle inlet temperature = 3001 K. 

Nozzle inlet velocity =5759 m/s. 

Nozzle inlet area = 0.0508 m. 

Design exit pressure equal to the ambient pressure at an altitude of 43.9 Km. 

Fig. 5.1 shows the nozzle contour obtained from MOC and Ref [4] respectively. Here, the 
contour obtained from MOC technique is slightly under estimated compared to the Ref 
[4] result for frozen flow. However, the trend is same for both nozzle contours. This 
slight variation in results is because of species consideration. In the present work, the 
nozzle inlet has 79 percent nitrogen and 21 percent oxygen. Gases other than these two 
are neglected due to their negligible amount. In Ref [4], the results are obtained by 
solving all governing equations of flow by iterative procedure and the gas has nitrogen, 
oxygen, water, hydrogen and hydroxide. 
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In hypersonic vehicle, wall velocity distribution is an important parameter. The velocity 
dish ibution is a function of static pressure at the inlet of the nozzle. It should be 
calculated more accurately to predict the value of both the force and moment Co- 
efficient Fig 5.2 shows wall velocity distnbution of both MOC and Ref [4] results. Here, 
both MOC and Euler s code results are almost same. Ninety percent of the velocity 
increment is gained in the initial portion of the nozzle. This is because of the sharp 
expansion comer placed at the inlet of the nozzle. 

The effect of design exit pressure or altitude on nozzle contour is given in Fig. 5.3. Here 
for non-dimensional pressure value of 0.5, the exit area is increased to very high value, 
because that coiresponds to very high altitude, where the exit Mach number is very high. 
In lower altitude, both the MOC and Ref [4] values are same. In this design, atmospheric 
conditions (pressure, temperature, density) are calculated using the standard formula 
given by NASA [27]. This may not be useful for very high altitude, where all the 
thennodynamic variables have lesser values. Therefore, this deviation in trend in the 
initial portion of the graph is mainly due to incorrect calculation of atmospheric pressure 
at very high altitude. The pressure at an altitude of 55 Eon is equal to 48.022 N/m^. 

Fig. 5.4 shows the drop in static pressure along the length of the nozzle. Here, the static 
pressure is constant for a while after the comer expansion fan, since there is no 
characteristic wave in that zone where the flow is uniform. Therefore, the nozzle contour 
IS a straight line up to the point where the first reflected characteristic wave meets the 
contour. However, in Ref [4], the pressure variation is smooth along the length of the 
nozzle. These results obtained from MOC are also verified with Euler’s code. 

5.2 Thermodynamic variables 

The thermodynamic properties such as pressure, temperature and density are calculated 
using the formulae given by National Aeronautics and Space Administration (NASA) for 
different altitudes. Fig 5.5, Fig 5.6 and Fig 5.7 shows variation of temperature, pressure 
and density with altitude. In Fig 5.5, the temperature plot consists of series of straight- 
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line segments and constant temperature segments. In troposphere, the temperature 
reduces as altitude increases because the radiation produced by the earth surface reduces 
as altitude increases. Hypersonic vehicles prefer altitude from upper stratosphere to very 
high altitude for flight. 

In Fig. 5.6, the pressure reduces exponentially as altitude increases. When altitude 
increases, mass of the substance in the atmosphere is reduced because of gravitational 
pull. For current generation flight, the materials available to withstand the load is in the 
range of 23,000 N/m — 95,000 N/m . In troposphere, the density of the air is very high 
compared to stratosphere. If Mach number exceeds five, the loads acting on the vehicle is 
very high compared to available maximum allowable load. Therefore, hypersonic flight is 
not possible in troposphere. 

In Fig. 5.7, density reduces exponentially with increase in altitude. Hypersonic flight is 
possible in stratosphere and above, because the density at this altitude is very less. 
However, hypersonic flight needs higher mass flow at inlet to sustain hypersonic flight. 
This is achieved by choosing larger stream tube area at the inlet with a series of oblique 
shock where the pressure loss is very less. Similarly, to expand the flow to hypersonic 
speed higher area at the nozzle exit is required. All these requirements are satisfied in 
hypersonic vehicle, when choosing the entire bottom surface of the vehicle as a thrust 
producing device. 

5.3 Parameters that affect the nozzle performance 

Fig. 5.8 shows the schematic diagram of single expansion ramp nozzle. Here, the uniform 
flow region is a zone where the flow properties are constant and has no waves Simple 
region has characteristic waves, which do not cross each other. If charactenstic waves 
cross each other, the region is called non-simple region. Fig. 5.9 to 5.13 shows the effect 
of different parameters that affect the nozzle performance for inlet Mach number 1 .5, 2.0, 
2.5 at a constant exit Mach number 5. Here the specific heat ratio value is 1 .24. 
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5.3.1 Area ratio 


Fig. 5 9 shows the effects of area ratios on nozzle performance. The area ratio is a 
function of exit Mach number. Here, the points 1, 2, 3 denote the area ratio required for 
exit Mach number 5 for inlet Mach number 2.5, 2.0, 1.5 respectively. For inlet Mach 
number 1 .5, the point after the initial expansion occurring near the wall (uniform region) 
has a Mach number of 2.84, with a corresponding area ratio of 3.5. This is a reasonable 
exit area with lesser length. The area ratio for exit Mach number 5 is 63.4, since this 
value is too large, hypersonic flight is not possible. In order to sustain hypersonic flight, 
the nozzle is truncated to get reasonable length and area. 

When the inlet Mach number increases to high values, the uniform flow region length 
behind the sharp expansion comer also increases. This leads to a straight-line shape in the 
initial portion of the nozzle with large flow turning angle. Therefore, for a very high inlet 
Mach number the nozzle becomes straight-line with large angle in the initial portion of 
the nozzle. 

5.3.2 Uniform flow region length 

Fig. 5.10 shows the effect of uniform flow region lengths on overall nozzle length of the 
vehicle. This uniform flow region length is a function of exit Mach number of the nozzle. 
If the exit Mach number increases to hypersonic Mach number range, then the unifomi 
flow region length is around 80 to 90 percent of the overall length of the nozzle, 
irrespective of the nozzle inlet Mach number. Slight increase in length causes both 
increase in weight and drag in hypersonic vehicle. Similarly if the nozzle inlet Mach 
number is very large (say M=5) then flap length also increases to a considerably large 
value. This produces additional drag. Therefore, in hypersonic flight both inlet and exit 
Mach numbers are important parameters. 
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5.3.3. Total flow-turning angle 


Fig. 5 1 1 shows the effect of total flow tuming angles on nozzle performance. Here, the 
point 3 denotes the total flow tuming required for exit Mach number 5 where the inlet 
Mach number is 1.5. From the graph, the total deflection angle required for Mach 5 is 87 
degree The flow tuming at the initial simple expansion fan comer is one 7 half of the total 
flow-tuming angle. Therefore, the initial expansion angle is 43.5 degree. This leads to a 
very large value of height, which produces a large drag. If the nozzle inlet Mach number 
is of higher value, (say M=2.5) then the initial expansion is around 26 degree. Therefore, 
this produces less drag compare to the previous case. In the case of higher nozzle exit 
Mach number requirements, the Mach number is increased to high value before it comes 
to nozzle inlet. 

5.3.4. Nozzle static pressure drop 

Fig. 5.12 shows the drop in static pressures along the length of the nozzle. Here, the point 
1 shows the drop in static pressure value at uniform flow region. 90% of the pressure 
energy is converted into kinetic energy at this zone. Therefore, the remaining 10% of the 
available pressure takes the entire remaining length of the nozzle for expansion. This 
drop in pressure is also a function of back pressure (ambient pressure). If the back 
pressure is high enough then the drop in pressure' in initial, expansion fan is reduced. 
Therefore, the amount of pressure energy converted to kinetic energy is also a function of 
exit Mach number. 

5.3.5. Stream thrust function 

Stream thrust function at a location is a measure of how much pressure energy is 
converted in to kinetic energy in the flow. So this function is used to measure the 
performance of the nozzle. Fig. 5.13 shows the stream thrast function with area ratio 
along the nozzle length. Here, point 1 denotes the stream thrust function at uniform zone. 
From this graph, more than 50% of the pressure energy is converted into kinetic energy in 
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the initial portion of the nozzle inlet because of sharp expansion comer. 60 percent of the 
stream thrust function value gives reasonable length. Therefore, for any supersonic 
nozzle design all the above parameters are to be considered to optimize the nozzle. 

5.4 Effect of Gamma on Nozzle exit parameters 

Specific heat latio is a sensitive parameter in hypersonic speed range. This value depends 
on mass concentration of each species, temperature, and pressure. The nozzle flow in 
hypersonic flight has very high temperature. At this temperature, ' either gas or gas 
mixture will begin to dissociate. Oxygen would dissociate at the temperature above 2200 
K. Nitrogen would dissociate at the temperature above 4500 K. Hence, up to 2500 K the 
gas is assumed as thermally perfect gas and the flow is frozen (the chemical composition 
of the exhaust is constant) with no dissociation. In this nozzle design, the chemical 
composition at the nozzle inlet is frozen. 

Specific heat capacity (Cp) is one of the sensitive parameters in hypersonic nozzle design 
In ordinary wind tunnel supersonic nozzle design, the specific heat ratio is assumed 
constant because in wind tunnels nozzle the maximum temperature is below 800 K. For 
temperatures up to 800 K, the gas is assumed calorifically perfect. Once the temperature 
increases above 800 K, the vibrational degree of freedom of the diatomic gas increases. 
This will reduce the value of specific heat ratio. 

Table 1 show the effect of specific heat ratio on nozzle parameters, where the input 
parameters are constant for all cases. The gamma value of 1.24 shows that the Mach 
number at exit is reasonably lesser compared to other two cases of gamma 1.33 and 1.4 
respectively. If gamma increases, all the exit parameters Mach number, exit area and 
length are also increases. The assumption of higher value of gamma leads to incorrect 
value of exit Mach number that leads to over estimated nozzle contour. If the assumed 
gamma value is less, then this leads to under estimation of the nozzle contour. Therefore, 


45 



the only way to predict the correct nozzle contour is variable gamma design, where 
gamma is function of static temperature throughout the nozzle length. 

The static temperature is not a constant value along the nozzle length. This reduces along 
the nozzle length when the flow begins to expand. The specific heat ratio is a function of 
temperature when the temperature exceeds 800 K. Therefore, the gamma value increases 
from a lesser value to a higher value along length of the nozzle because of reduction in 
vibrational degree of freedom, which makes the constant gamma assumption invalid for 
hypersonic nozzle design. The variable gamma is calculated along the length of the 
nozzle using the formula given in Ref [4]. Apart form this variable gamma of gas, the 
fuel used in the combustion chamber is also taken into account. Here, the combustion 
takes place in the combustor is complete. So, the stoichiometric equation is used to 
calculate the exact composition of the exhaust. The specific heat capacity value is 
calculated from the formula given in Ref [5]. The gamma value for propane fuel is also 
calculated from this formula and the corresponding results are given in the table. 

Fig. 5.14 shows the effect of gamma value on nozzle contour for gas and propane fuel 
with both constant and variable specific heat ratio. Here the design exit pressure is equal 
to the ambient pressure at an altitude of 30000 m. The other input parameters are given 
below. 

Nozzle inlet Mach No = 1.500 
Nozzle inlet static temperature = 1800.000 K 
Nozzle inlet static pressure = 306551.6401 nW 
The inlet height or area of nozzle = 0.600 m 


Here, in the constant gamma case the gamma value of 1.33 gives good approximate value 
for exit area compared to other two cases. The nozzle is truncated to a lesser length where 
this approximate gamma gives best contour. In variable gamma gas, the air contains 79 
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percent nitrogen and 21 percent of oxygen. Since the gas composition is known, the 
specific heat ratio value for air at different temperature is calculated. Similarly, the 
piopane fuel has three carbon atoms and eight hydrogen atoms. This fuel reacts with air 
and gives the products. This product composition is calculated using stoichiometric 
equation. Here, the exit area of propane fuel is slightly higher compare to variable gamma 
of gas Thei eforc, the correct value of exit area is predicted when the fuel is also taken 
into account. 


5.5 Effect of back pressure 


The exit area of the nozzle is a function of both nozzle inlet static pressure and ambient 
pressure. Fig. 5.15 shows the effect of back pressure or altitude on exit area or nozzle 
contour of the nozzle. The nozzle inlet pressure is constant for all three cases. The exit 
area for constant gamma value is large compared to variable gamma at very low altitude 
or high-pressure region. However, the difference in exit area for propane and variable 
gamma of gas arc very less. For low-pressure region or very high altitude (>50 km) all 
three have considerable difference in exit area. Therefore, for hypersonic flight speed the 
specific heat ratio is an important factor. 

If the back pressure increases then the Mach number at the nozzle exit reduces because of 
high backward force. This reduction in Mach number causes the reduction in exit area. 
For con'cctly expanded flow at that, altitude needs less exit area. The entire nozzle 
contour is modified when the back pressure increases. The higher exit area leads to over 
expansion and lesser exit area causes under expansion. In both the process, loss of thrust 
takes place, therefore the given nozzle contour is suitable only for a particular altitude. In 
higher back pressure region, the variation of exit area with pressure is less compared to 
lesser back pressure region (high altitude). The variation of nozzle contour is large in 
high altitude where the hypersonic flight is possible. The cruising altitude must be 
constant for hypersonic vehicles. Therefore selecting a cruising altitude for hypersonic 
flight is a difficult task for designers where all the parameters must be considered. 
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5.6 Variation of Nozzle length with number of characteristics 


The supersonic flow is wave-dominated flow with infinite number of characteristics 
moving both right and left sides along the flow. However, it is not possible to consider 
infinite number of characteristics in computation. The number of characteristics is 
maximized depends on computer s efficiency and speed. Therefore, for simplicity finite 
number of characteristics is considered. The accuracy of the problem depends on number 
of characteristics. 

Table 2 shows the variation of total length of the nozzle with number of characteristics 
for constant gamma for given inputs. If the number of characteristics is less then the 
nozzle length is also less. If the number of characteristics increases then the total length 
of the nozzle is increased. However, after certain number of characteristics the length is 
almost remains constant. Therefore, the less angle (<=1 degree) between the 
characteristics gives results with good accuracy. For lesser angles, impact on nozzle 
length is very less. Table 3 shows the variation of nozzle length with number of 
characteristics for variable gamma. 

The reduction in number of characteristics leads to increase in the hngle between two 
successive Mach lines. This large angle shows shock wave instead of Mach wave because 
the changes in properties across the Mach wave are very small. The reduction in number 
of characteristics leads to shock instead of Mach line in the flow. Therefore, the number 
of characteristics is assumed as large as possible in the total deflection angle. This leads 
to accurate solution. 

5.7 Results for test case 

The following results are for two-dimensional nozzles. In Method of characteristic 
technique both two dimensional and axi-symmetric nozzle are s^e because each 
characteristic passes and reflects on its own plane. Therefore, the only difference in axi- 
symmelric flow the input is radius, whereas in 2-d nozzles it is the inlet area. Here, the 
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oulpul file from the program is given where inlet Mach number, inlet area, inlet pressure 
and temperature are inputs along with altitude. The corresponding results in graphical 
form arc also given figures from Fig. 5.16 to Fig. 5.31. 

The one general output file for variable gamma for test case is given below: 

Atmospheric conditions for 30000.000 m altitude 
Atmospheric tcmperature= 23 1 .490 K 
Atmospheric pressure - 1 161.181 N/m^ 

Atmospheric density = 0.017 kg/m^ 

INPUTS 

Nozzle inlet Mach no = 1.500 
Nozzle inlet Mach angle = 0.7297 rad 
Nozzle inlet static temperature = 1800.000 K 
Nozzle inlet total temperature = 2410.243 K 
Nozzle inlet static pressure = 306551.640 N/m^ 

Nozzle inlet total pressure = 1081483.876 nW 
The inlet height or area of nozzle = 0.600 m 
Separation angle in radians = 0.0017 rad 

OUTPUTS 

Nozzle inlet local velocity of sound = 821.753 m/s 

Nozzle inlet velocity = 1232.629 m/s 

Required Mach no at exit = 5.485 

Exit Mach angle = 0.183 rad 

Exit temperature = 346.089 K 

Exit velocity of sound = 373.297 m/s 
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Exit velocity ^ 2047.732 m/s 
Throat area ~ 0.504 
Area ratio a/a* at inlet=1.189 
Area ratio a/a* at exit =37.580 
Exit area of the nozzle =31.598 
Total flow turning angle = 1.201 rad 
Total How turning angle = 68.852° 

Total number of characteristics = 345 

First characteristics angle = 0.000462 rad 

Flow turning angle at inlet = 0.000 rad 

Mass flow rate 452.670 Kg/s 

Nozzle overall length = 1 1 8.771 m 

Unifonn flow region length = 102.246 m 

Flap length required for desired exit Mach number = 16.524 m 

Fig. 5.15 shows the nozzle contour for inlet Mach number 1.5 where the exit Mach 
number is 5.5. From this, the inference is that, the initial expansion sharp comer angle is 
very large because the exit Mach number is very large also; the overall length of the 
nozzle is equal to 118m. Therefore, from the above data, the flap length required for very 
high Mach number is very small compared to the overall length. The initial exit area at 
the inlet for this nozzle is very large compare to large inlet Mach number conditions. This 
large area at the initial portion causes large amount of drag. During such situations, 
conical nozzle yields good results compare to contour nozzle. In order to get such a high 
Mach number at exit through contour nozzle design, the nozzle inlet Mach number also 
should be high. This can be achieved by increasing the nozzle inlet area. Finally, the 
length of the nozzle is optimized using stream thrust function. 
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Fig 5.16 and Fig. 5.17 show the actual and enlarged view of intermediate points. This 
gives an idea of how the flow properties vary across the charactenstic wave The flow 
piopeitics vaiiation is large in the initial portion of the nozzle compared to aft region of 
the noz/lo This means the change in properties or expansion takes rapidly in the initial 
poition ol the nozzle. Heie the characteristic wave is straight line between any two 
points. However, in reality this is not a straight line because the local temperature is 
varying at every point in the flow. 

Fig. 5.19 and Fig. 5.20 show the drop in static pressure and static temperature along the 
length of the nozzle respectively. Here the drop in pressure and temperature is very large 
in the initial portion compared to aft region of the nozzle. Here the figures show that at 
the nozzle inlet the drop in pressure is very large because of sharp expansion comer used 
at the inlet of the nozzle. These changes in properties are given in Fig. 5.17 also. 
Therefore, wc can truncate the nozzle at certain location so that the reasonable exit Mach 
number is obtained at the exit. This can be optimized using stream thrust function. 

Fig. 5.21 and Fig. 5.22 show the variation of stream thrust function with area ratio and 
length of the nozzle respectively. Here most of the available pressure energy is converted 
into kinetic energy in the initial portion of the nozzle. In hypersonic flight drag is a main 
parameter, which depends on both exit area and length of the nozzle. Therefore, both 
length and height are optimized using this parameter. Five times the inlet area of exit area 
gives reasonable length and height. The length of the nozzle also depends on the exit 
Mach number requirement. 

Fig. 5.23 and Fig. 5.24 show the axial and normal force coefficient respectively. Here the 
axial and the nomial force coefficients are nomialized using exit dynamic pressure. Both 
axial and normal force coefficients are very large at the initial portion of the nozzle and 
this variation is small as the length increases. This axial force leads to additional lift in 
the nozzle surface, which produces pitching moment. • 



Fig. 5.25 and fig. 5.26 show the axial and normal moment coefficient along the length of 
the noz/Jc respectively. The stability parameters depend on the value of both axial and 
nonnal moment coefficients. The incorrect prediction of these coefficients leads to 
unstable configuration. This predicted value of both axial and normal force coefficients is 
suitable foi a particular altitude. Therefore, in order to fix certain static margin for 
stability, these values must be calculated precisely. 

Fig. 5.27 to Fig. 5.29 shows the variation of different types of thrust along the nozzle 
length. F'roni h'ig. 5.28, at the nozzle inlet the pressure thrust is large because at the inlet 
available pressure energy is large. This pressure energy is converted into kinetic energy 
when the flow expands down stream of the flow. Fig. 5.29 shows the variation of net 
static thrust along the length. The net thrust value at the inlet is very less. However, after 
the shapp expansion this thrust is increased to a larger value. After that, the increase in net 
thrust is very less. Therefore, the length is truncated after the sharp expansion depending 
on required Mach number at the exit. 


Fig. 5.30 and Fig. 5.31 shows the truncated nozzle contours for given inlet conditions. 
Here the nozzle is truncated where the stream thrust function is 70 percent at the nozzle 
wall. The overall nozzle length for truncated nozzle is 3.8147 m and the exit area is 
2.0945 m. The flap length required for tmneated nozzle is 2.7 m. This gives the 
reasonable length and exit area compared to the full length and area of the nozzle. The 
exit Mach number for truncated Nozzle is 3.3818. 

In the case of full nozzle, the exit Mach number is 5.5 where the length of the nozzle is 
118.771 m and area is 31.599 ml This higher value of both length and area is not 
possible in hypersonic vehicles. Therefore, the nozzle is truncated depending on the 
required exit Mach number of the vehicle. The stream thrust function value around 60 
percent, gives the reasonable exit area and exit Mach number. However, the truncated 
nozzle has flow with some angularity. 
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5.8 The results obtained from Euler’s code 


The icsulls ohtiUned Tioiti the NtOC ^re verified with Euler^s code results of inviscid 
now In Euler’s code, the contour of the Scramjet nozzle is designed using the contour 
obtained from present MOC method, and then the flow through the nozzle is analyzed by 
solving all the governing equations for given input parameters. Here, the specific heat 
ratio gamma is constant throughout the length of the nozzle. The MOC and the Euler’s 
code results at exit are shown in Table 4. 


The input and the output from the Method of characteristic numerical method for 
constant gamma are given below. 

Free stream conditions 

Altitude =30000 m 

Free stream temperature =23 1 .49 K 

2 

Free stream pressure =1161.180 N/m 

Free stream density =0.0 1 74 Kg/m^ 


Inlet Conditions 

Nozzle inlet Mach no 
Inlet temperature 
Inlet pressure 
Inlet density 
Inlet velocity 
The inlet area 


=1.5 

=1800 K 

=306551.640 NW 
=0.593 Kg/m^ 
=1243.352 m/s 
= 0.600 m 
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Exit Conditions 


Nozzle exit Mach no 
Exit temperature 
Exit pressure 
Exit velocity 
Exit area of the nozzle 


= 5.205 
= 451.256 K 
=1 161.180 NW 
= 2160.1158 m/s 
=22.857 m^ 


Fig. 5.32 and Fig. 5.33 show the Mach number and pressure distribution inside the nozzle 
for hypersonic vehicle. Fig. 5.32 gives a good Mach number distnbution compared to 
Method of characteristic results. Here each color zone corresponds to one particular 
Mach number. These zones are similar to the characteristics net explained in Fig. 5.18 
where the Mach wave is straight line in between two points. However, in reality it is not a 
straight line. For simplicity, the Mach waves are assumed straight line in between two 
intermediate points. This Mach number distribution is for 5000 iterations and it takes 1 
full day for 5000 iterations. However, this can be solved with in few minutes using 
Method of characteristic technique. 

The drop in static pressure is given in Fig. 5.33. Here the drop in pressure is large in the 
initial portion of the nozzle. The aft portion of the blue color region shows the uniform 
flow in the nozzle where the pressure is constant. Fig. 5.34 shows the comparison of 
Mach number distribution along the nozzle wall for both MOC and Euler s code results. 
Here the MOC result shows the gradual variation of Mach number along the wall 
compare to Euler’s code computational results. 

In the initial portion, both curves are same because at this portion the expansion takes 
place rapidly. Therefore, the difference in value at each grid is large. However, in the 
down stream the expansion is small compared to initial portion. Large number of grids in 
the aft portion leads to waste of time. Less number of grids in the aft portion of the nozzle 
gives uniform variation of Mach number and flow properties. 
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5.9 Important fact 


The designed nozzle is single sided and single expansion ramp nozzle. Here a flap is used 
for reflecting the characteristics and it is used to vary the exit conditions for over and 
under expanded cases. In the case of double-sided symmetric nozzle one-half of the 
nozzle inlet area in taken in to consideration for contour design. Here, the flow is 
symmetric in both sides along the length of the nozzle. This is not used m hypersonic 
vehicle because this produces enormous drag compare to net thrust. Therefore, in order to 
design a symmetrical nozzle the inlet area of the program should be one-half of the total 
inlet area of the nozzle. 
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Fig. 5.35 NASP (USA) hypersonic vehicle 

The current generation vehicles are mostly single sided conical nozzle. The conical 
nozzle exit has the flow with some angularity. Therefore, considerable loss in thrust takes 
place in such nozzles Fig. 5.35 shows the conventional National Aero Space Plane 
(NASP) hypersonic vehicle, which is a single sided expansion nozzle. 
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The contour nozzle has considerable axial flow at the exit of the nozzle because the 
nozzle is truncated at certain location. The flap is used for reflecting the characteristics or 
Mach waves. During off design conditions the flap can be moved up and down to adjust 
the flow in the nozzle exit to a correctly expanded flow. However, deflecting the flap at 
hypersonic Mach number is a very difficult task. The following Fig 3.35 shows the single 
expansion ramp nozzle. 



Fig. 5.36 Single expansion ramp nozzle 

Here at the nozzle exit the flow is under expanded because the initial correctly expanded 
flow nozzle is truncated at some reasonable length. Therefore, it should have some 
expansion waves at the exit. However, change in property across this wave is small 
compared to initial expansion. It too has flow with certain angularity at the exit. 
However, this nozzle is efficient compared to conical nozzle because at the initial portion 
contour nozzle has more area for expansion whereas in conical nozzle, this area is less. 
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CHAPTER 6 


CONCLUSION AND REMARKS 

The method of characteristics computer program for designing and analyzing the 
Scramjet nozzles has been developed. This program has also incorporated chemical 
composition of air and fuel used in the combustion chamber for calculating the specific 
heat ratio, where complete combustion is assumed. The calculated chemical composition 
in the nozzle inlet is assumed constant throughout the length of the nozzle. The 
temperature and pressure along the wall at different locations are calculated, using which 
the force coefficients are calculated. Provision is provided for calculating the nozzle 
contour at different altitudes. This code can also be used for axi-symmetric nozzles. The 
code was used here to design a generic two-dimensional hypersonic nozzle based on 
arbitrarily selected flight conditions. This program is used for static temperature range up 
to 2500 K with best accuracy. However, it can be used for temperature up to 3000 K with 
a slight lesser accuracy. An evaluation of the various effects studied leads to the 
following conclusions. 

1. Scramjet nozzles have good performance estimates when gamma has variable values 
rather than having various assumed constant values. Here flow is assumed to be frozen 
from the nozzle inlet to exit. The exhaust gas used in the combustion chamber also 
influences the performance of the nozzle. 

2. Here, the chemical composition is assumed constant throughout the nozzle. If the 
temperature is very high, the chemical composition varies across each wave while 
expansion. Therefore, up to 2500 K, the gas is assumed as a function of temperature. 

3. In the case of Scramjet nozzles, the separation angle for each characteristic should be 
less than one degree. In our case, 0.25° gives good results. 
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4 If the Mach number is of low supersonic at the inlet, while at the exit the Mach number 
IS of hypersonic, then the initial portion of the nozzle has a larger height. In such 
situation, two single expansion ramp nozzles are preferred. 

5. Most of the available pressure and temperature are converted into kinetic energy in the 
initial portion of the nozzle. Therefore, stream thrust function value of sixty percent gives 
reasonable height or area and length for hypersonic vehicles. 

6. The angularity of flow in the nozzle wall region produces additional lift in the nozzle 
portion, which produces pitching moment about the centre of gravity of the vehicle. 
Therefore, these parameters are used to calculate vehicle static margin or trim 
requirements of the vehicle. 

7. Momentum and pressure thrust along the axial location of the nozzle at different 
locations are calculated. The net thrust is calculated by adding up the momentum and 
pressure thrust. 
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CHAPTER 7 


FUTURE SCOPE OF WORK 

The given Method of characteristics nozzle contour design is suitable for correctly 
expanded flow, but practically this is not feasible The given design is for one particular 
cruising altitude only. Therefore, at all other altitudes this vehicle operates in off design 
conditions. The following are some ideas for future scope of work. 

1 The flow through the Scramjet nozzle is assumed as thermally perfect gas. However, 
the temperature after 2250 K the gas starts to dissociate. The chemical composition of gas 
across each characteristic wave differs when the flow expands along the length of the 
nozzle Therefore, the chemical composition across each characteristic must be calculated 
con-ectly to predict the performance parameters 

2. Combustion taking place in the combustion chamber of hypersonic vehicle is assumed 
complete. However, it never happens in reality. Here the small portions of all other 
compositions are also neglected. Therefore, the exact chemical composition must be 
calculated. 

3. This given Scramjet nozzle design is used for inviscid flow. However, in real nozzles, 
there should be a boundary layer on the wall. The boundary layer thickness can be 
calculated for viscous flow to predict the exact nozzle contour. This program can be 
modified for both chemically reacting flows with or without chemical non-equihbnum. 

4. This nozzle design is suitable for correctly expanded flows. However, except the 
particular altitude the flow is either over expanded or under expanded. The under 
expanded flow can be modified when deflecting the flap in the downward direction. The 
over expanded flow can be modified by deflecting the flap in the upward direction. 
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results 

TABLE 1 

Variation of Exit condition with Specific Heat Ratio 

INLET CONDITIONS 


Altitude of flight 30000 m 

Atmospheric T emperature 23 1 .49 K 

Atmospheric Pressure 1 1 6 1 . 1 8042 N/ra^ 

Atmospheric Density 0.017476 kg/m^ 

Inlet Mach number 1 .500 

Nozzle Inlet Static Temperature 1 800.00 K 

Nozzle Inlet Static Pressure 306551 .640140 N/m 

Inlet Area of the Nozzle 0.60 m^ 


2 


Specific heat ratio 

(Gamma) 

Nozzle Outlet Conditions 

Gamma=l .4 

Exit Mach number 5.537470 

Exit Temperature 365.919597 K 

Exit Velocity 2123 . 289318 m/s 

Exit Area of the Nozzle 19 . 34 6021 m^ 

Nozzle Overall Length 125.348746 m 

Required Flap Length 19.981207 m 
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Gamma=1.33 


Exit Mach number 5.204734 
Exit Temperature 451.256437K 
Exit Velocity 2334.298584 m/s 
Exit Area of the Nozzle 22. 857215 m^ 
Nozzle Overall Length 139.261340m 
Required Flap Length 22.511033 m 


Gamma=1.24 


Exit Mach number 4.775643 
Exit Temperature 611.751399 K 
Exit Velocity 2228.287071 m/s 
Exit Area of the Nozzle 29 . 004528 m^ 
Nozzle Overall Length 162.186408 m 
Required Flap Length 26.740636 m 


Correct prediction 
Variable gamma 


Correct prediction 
Variable gamma 
for propane 


Exit Mach number 5.485522 
Exit Temperature 346.089108 K 
Exit Velocity 2 04 7. 7 3 2415 m/s 
Exit Area of the Nozzle 18.959363 m^ 
Nozzle Overall Length 118.771246 m 
Required Flap Length 16.524270 m 

Exit Mach number 5.321639 
Exit Temperature 379.071917 K 
Exit Velocity 2072. 658005 m/s 
Exit Area of the Nozzle 20.431952 rci 
Nozzle Overall Length 123.218631 m 
Required Flap Length 16.540835 m 
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TABLE 2 


Variation of nozzle length with number of characteristics for constant 

gamma 


INLET CONDITIONS 


Altitude of flight 
Atmospheric Temperature 
Atmospheric Pressure 
Atmospheric Density 
Inlet Mach number 
Specific Heat Ratio 
Nozzle Inlet Static Temperature 
Nozzle Inlet Static Pressure 
Inlet Area of the Nozzle 


30000 m 
231.49 K 
1161.18042 N/m^ 
0.017476 kg/m^ 

1.500 

1.3300 kJ/kg 
1800.00 K 

306551.640140 nW 
0.60m^ 


EXIT CONDITIONS 

The exit Mach number 
Exit Temperature 
Exit Velocity 
Exit area of the nozzle 


5.204734 
451.256437 K 
2160.115896 m/s 
22.857215 m^ 
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Number of Characteristics 

Exit conditions 

Total number of characteristics 20 

Nozzle overall Length 132.666028 m 

Flap Length 15.915722 m 

Total number of characteristics 39 

Nozzle overall Length 135.729271 m 

Flap Length 18.978964 m 

Total number of characteristics 76 

Nozzle overall Length 137.584508 m 

Flap Length 20.834202 m 

Total number of characteristics 151 

Nozzle overall Length 138.610820 m 

Flap Length 21.860513 m 

Total number of characteristics 375 

Nozzle overall Length 139.261340 m 

Flap Lengtli 22.51 1033 m 
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TABLE 3 


Vai iation of nozzle length with number of characteristics for variable 

gamma 

INLET CONDITIONS 


Altitude of Flight 

30000 m 

Atmospheric Temperature 

231 49 K 

Atmospheric Pressure 

1161.18042 N/m^ 

Atmospheric Density 

0.017476 kg/m^ 

Inlet Mach number 

1.500 

Specific heat ratio 

1.33KJ/kg-K 

Nozzle Inlet Static Temperature 

1800.00 K 

Nozzle Inlet Static Pressure 

306551.640140 nW 

Inlet Area of the Nozzle 

60.000 cm^ 

EXIT CONDITIONS 


The exit Mach number 

5.204734 

Exit Temperature 

451.256437 K 

Exit Velocity 

2160.115896 m/s 

Exit area of the nozzle 

22.857215 m^ 
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Number of Characteristics 

Exit conditions 

Total number of characteristics 19 

Nozzle overall Length 113.082561 m 

Flap Length 11.481438m 

Total number of characteristics 36 

Nozzle overall Length 1 15.830752 m 

Flap Length 13.788870 m 

Total number of characteristics 70 

Nozzle overall Length 1 17.390839 m 

Flap Length 15.218479 m 

Total number of characteristics 139 

Nozzle overall Length 118.238067 m 

Flap Length 16.016498 m 

Total number of characteristics 346 

Nozzle overall Length 118.771246 m 

Flap Length 16.524270 m 
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TABLE 4 


COMPARISON OF RESULTS 


PARAMETERS 

MOC RESULT 

EULER’S CODE 
RESULT 

PERCENTAGE OF 
ERROR 

Mach number 

At exit 

5.205 

5.189 

1.067 

Static pressure at 
exit 

1161.168 N/m2 

892.245 N/m^ 

.0877 
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Non dimensional 
velocity(v/v_e) 


Fig. 5.1 Nozzle contour for frozen flow 



Fig. 5.2 Wall Velocity profile 





Fig. 5.3 Effect of Exit Pressure 



Fig. 5.4 Drop in Static pressure along the wall 
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Zone VI Length to Height 



0 2 4 6 8 

Nozzle Exit Mach no 


Fig. 5.9 Effect area ratio on nozzle performance 



Fig. 5.10 Effect of uniform flow region length on nozzle length 
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Fig. 5.11 Effect of Total flow turning angle on nozzle performance 



0 2 4 6 8 

Nozzle Exit mach number 


Fig. 5.12 Effect of static pressure on nozzle length 
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Fig. 5.13 Effect of Stream thrust function along the nozzle height 
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Effect of Back pressure 



Fig. 5.15 Effect of back pressure on exit area 



!■ 




Fig. 5.17 Intermediate points 




Non dimensional Non dimensional 

pressure(p/pj) height(h/h_i) 


Intermediate points 
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Non dimensional length(l/hj) 


Fig. 5.18 Intermediate points Enlarged view 


Pressure drop 



0 50 100 150 200 250 


Non dimensional length(l/h_i) 
Fig. 5.19 Pressure drop 
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Fig. 5.23 Axial force coefficient 
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Fig. 5.26 Normal Moment coefficient 
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Fig. 5.28 Pressure 


thrust along the length 
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110 1 
Non Dimensional Length (l/hj) 


Fig. 5.30 Truncated nozzle contour 


Nozzle Contour 






Results derived from Euler code 



Fig. 5.32 Variation of Mach number along the length 
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Fig. 5.33 Variation of pressure along the length 
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Fig. 5.34 Comparison of Mach number along the wall 
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